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Probing the birth of fast rotating magnetars through high-energy neutrinos 
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We investigate the high-energy neutrino emission expected from newly born magnetars surrounded 
by their stellar ejecta. Protons might be accelerated up to 0.1-100 EeV energies possibly by, e.g., the 
wave dissipation in the winds, leading to hadronic interactions in the stellar ejecta. The resulting 
PeV-EeV neutrinos can be detected by IceCube/KMSNet with a typical peak time scale of a few days 
after the birth of magnetars, making the characteristic soft-hard-soft behavior. Detections would 
be important as a clue to the formation mechanism of magnetars, although there are ambiguities 
coming from uncertainties of several parameters such as velocity of the ejecta. Non-detections would 
also lead to useful constraints on the scenario. 



PACS numbers: 95.85.Ry, 98.70.Sa 



I. INTRODUCTION 

Magnetars are neutron stars endowed with the 
strongest magnetic fields known in the universe, B ~ 
;^Qi4-i5 Q ^fQj. reviews see [i|). It is believed that - 10 % 
of young galactic neutron stars possess such strong fields 
(for reviews see [l[ ) . Although the precise origin of these 
strong fields is uncertain. It has been argued that their 
amplification occurs via the dynamo mechanism during 
the Kelvin-Helmholtz cooling time, ^kh ~ 10 — 100 s, in 
the protoneutron star (PNS) phase 0. The dynamo ef- 
ficiency is partly determined by the initial rotation rate 
of rii = 2'K/Pi, and analytical estimates suggest that the 
formation of global strong magnetic fields might require 
Pi ~ 1 ms at birth. This will affect the dynamics of the 
PNS wind, by providing a significant reservoir of rota- 
tional energy comparable to that of the accompanying 
supernova (SN) explosion. The strong magnetic fields 
can facilitate the SN explosion, and may be responsible 
for the more powerful sub-class of SNe known as hyper- 
novae (HNe), related to long gamma-ray bursts (GRBs) 
(e.g. [3] for reviews), which together with the PNS wind 
acts as a piston on the compressed stellar ejecta, i.e. the 
young expanding SN remnant (SNR) (3]. 

Newly born magnetars with rapid rotation rates may 
be efiicient ultra-high-energy cosmic ray (UHECR) accel- 
erators 0, 0] . Particle acceleration in neutron stars has 
been considered at both the polar or outer gap inside the 
magnetosphere, and near the wind zone [U, Q • If ions are 
accelerated as well as electrons, up to sufficiently high en- 
ergies, they can produce high-energy neutrinos via, e.g., 
the jry or pp reaction. Based on the polar gap models, 
Rcfs. discussed neutrino production via the P7 re- 
action between ions and surface x rays, or by curvature 
pion radiation of ions. In this work, we investigate the 
high-energy magnetar neutrino emission resulting from a 
different scenario sueg estedbyRef. Q , where cosmic-ray 
acceleration is attributed to the wake-field acceleration 
mechanism beyond the light cylinder. We consider the 



pp and p7 interactions of cosmic-ray ions with cold SNR 
nucleons and thermal photons. In the case of fast ro- 
tating magnetars, we find that the resulting neutrino en- 
ergy fluence peaks around t a, few days due to hadronic 
cooling of mesons and muons. Such high-energy neutri- 
nos can be detected by future km^ telescopes such as 
IceCube/KMSNet [| in a few years or if magnetars are 
born at < 10 Mpc. 



II. THE MODEL 

In the dynamo scenario, newly born magnetars have a 
large rotational energy 0. During <kh the PNS winds 
would be thermally neutrino-driven or magnetically dom- 
inated but subrelativistic, and the spin down rate may 
be enhanced by neutrino-driven mass loss After ^kh, 
the winds become magnetically dominated and relativis- 
tic, similarly to the case of pulsar winds. The rota- 
tional energy is extracted by the Poynting flux and grav- 
itational waves Using the vacuum dipole formula, 
the rotational energy loss rate by the magnetic wind at 

t{> Tem ^ lO'-^ S /45 

A^33^ 4 ) ill the stellar frame can 
be estimated as L{t) ~ 6.1 x 10^^ erg s~^ I^^ ii'^^ 
where ^ = i^NS^NS - 0-5 x 10^^ G cm^ ^ns,15-Rns,6 
the magnetic dipole moment and the moment of inertia / 
is set to 10"*'^ g cm^ throughout this work. The particles 
are expected to gain a fraction of the magnetar rotational 
energy during their acceleration, by tapping a fraction of 
the open field line voltage on their ways from the mag- 
netar to the outside region. The maximum cosmic-ray 
energy accelerated at t(> Tem) is 



E^'{t) = ?/Ze$,„ag ^ 2.0 X 10^0 cV Z ry_i 145 ^^33' ^l'- 

(1) 

where 77 parameterizes the uncertainties in the utiliza- 
tion of the potential drop. In polar gap models, the 
parallel electric fields would be significantly screened in 
very young pulsars, implying 77 ^ 0.1 [H, Neverthe- 
less, as in pulsars, a significant fraction of the Poynt- 
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ing energy could be converted to the kinetic energy well 
outside the light cylinder, via mechanisms such as surf- 
riding acceleration. Ref. Q suggested UHECR accel- 
eration by the wake-field acceleration in the equatorial 
wind, where it was argued that, for an oblique rotator, 
much of the Poynting flux might be tied up in waves 
and the crinkled frozen-in current sheets dissipate around 
Tdiss ~ 10^^'^{c/n). If wave emission is the relevant 
dissipation process, it might form large amplitude elec- 
tromagnetic waves pushing ions by the ponderomotive 
force, Fpond ~ mc^{Z e5 B / mcVl) [1Q|. Hence, the work 
on a particle moving a distance / through the wave is 
^pond^ ~ ??^e<I>mag, as loug OS 5B ^ B and -q = l/r ~ 0.1. 
Following Ref. [a| , we hereafter assume that such prompt 
cosmic-ray acceleration mechanism is in operation, and 
use Eq. ([T]). The exact nature of the dissipation and ac- 
celeration mechanisms in the wind is currently uncertain 
and a detailed study is beyond the scope of this paper. 
Approximately, since we may expect that the ion injec- 
tion rate around the equatorial sector is the Goldreich- 
Julian rate [ll|, the cosmic-ray spectrum can be written 
as [i,[l2 



dN 
dE 



1 



(2) 



8 ZefiE{l + E/EG) 

where Eg = {b/72){riZe^^ / GPt^) and, for simphcity, we 
have assumed that all the particles accelerated at t have 
E^'{t). Note that since the acceleration is expected to 
occur promptly, energy losses can be neglected. Hence 
the adiabatic and radiation losses in the wind become 
irrelevant since accelerated ions will not be coupled to 
the fields and their curvature radius is large enough Q. 

III. THE NEUTRINO SPECTRUM AND FLUX 

A newly born magnetar will be surrounded by the 
young SNR, separated by a cavity evacuated by the 
wind or rapidly expanding SN shock. The particle ac- 
celeration in the wind occurs around race ~ ''diss ~ 
10^°"'' cm ^i^^t^J^ (hereafter we consider the cosmic-ray 
ions to be protons [HI). The wind termination shock 
and the SN shock radii are both larger than race for suffi- 
ciently late times as considered here. Thus, a significant 
fraction of cosmic rays will interact with the stellar ejecta, 
unless the latter is punctured or disrupted by the wind 
itself (see below for general cases). First, we consider 
the interaction between cosmic rays and SNR nucleons 
assuming that cosmic rays are emitted isotropically, to 
obtain conservative results. 

The cosmic rays will interact with the SNR via pp reac- 
tions, producing mesons. When the (magnetar-powered) 
SN shock has a high velocity of /3snc, the eff'ective op- 



tical depth for the pp reaction is /j 



pp ~ KppCTppJlp^^^ 



5.7 X 10^ A/sN,i/3sN 



-1^4 ' 



where k, 



pp 



0.5 



10^25 



0.6, app ~ 
(in the 100 PeV range) [H, and Asn 
''SN ~ l3sNct ~ lO^'^-^ cm /3sN,-ii4- The meson produc- 
tion efficiency is estimated as /mcs ~ niin(l,/pp). Since 
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FIG. 1: The {v/j, + i'fj,) fluence from a newly born magnetar at 
5 Mpc, at different time intervals. The fluence peaks around 
t ~ 2 days because hadronic cooling of mesons and muons is 
important at earlier times, while the amount of cosmic rays 
decreases with time. Thick/thin lines are for the cases with- 
out/with the radiation field. Vacuum neutrino oscillations are 
considered. 



the effects of magnetic fields in the SNR can typically be 
neglected, the interaction time is tint ~ Asn/c. 

The resulting charged mesons decay into neutrinos via 
TT^ — > -|- i^e(f'e) + i^fj. + v^. Thc ncutriuo spectrum 
roughly follows the cosmic-ray spectrum, but the high- 
energy spectrum is modified when mesons and muons 
cool down before they decay. In our cases, the in- 
elastic T^pj [ip collision can be relevant. When ]^p k, 
iint/^Trp ^ 1, the pion neutrino flux is suppressed as 
/sup ^ niin(l, (t7rp/77r'''7r)), where the break energy is 
Shad ^ 25 {t^plT^)m^c^ ^ 32 TeV M^^ ^^01,^ _^t\. 
'^njii is the proper life time of charged pions/muons and 
^ttp/mp ~ {K'-^p/^,pC^Tvp/ti.pnpc)~^ is their hadronic coohng 
Since the neutrino spectrum is proportional to 



time. 



/mes/sup, its rougli cxprcssion for dN / dE cx _B ^ is 



{E^/E^^-'^t 

(£;./£;had)i- 



(for < i^i-d) 
(for E^^'^^ < E^< i^;'"^'^) 

The resulting neutrinos can propagate in the SNR 
matter without significant attenuation because r^p w 



T^pUpAsN 



0.017 E°^i%MsN,if3, 



SN, 



-1*4 



« 1. 



We have performed detailed numerical calculations to 
evaluate neutrino spectra and fluxes through the method 
of Refs. [l3|. The calculations are performed during a 
time tint, taking into account the high meson multiplicity 
of the high-energy pp reaction, based on the SYBILL 
code fisj . and hadronic cooling of pious and muons with 
the approximated cross sections of a^^p ~ 5 x 10 



-26 cm2 

and cr^p ~ 2 x 10"'^** cm^ [l^. We neglect contributions 



from meson production via irp/fip processes which can 
affect the spectra by factors 0(1), since their influence 
is modest and only at relatively early times, t < day. 
As shown later, the detailed spectra will be somewhat 
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different from Eq. ([3]) especially at high energies, due to 
the high meson-multiplicity of the pp reaction and the 
accumulation of cooled mesons and muons. The only 
necessary input quantities are the cosmic-ray flux and 
target nucleon density. The former is given by Eqs. ([T]) 
and ([2]) and the latter by Msn and /3sn- 

In Figs. 1 and 2, the resulting spectra and light 
curves are shown for M^n = Msn.i = IOMq and 
/3sN = 0.1. We can see that the neutrino energy flu- 
ence becomes maximal around t ^ 2 days, because of 
hadronic cooling of pions and muons at earlier times 
and a decrease of E'^^{t) oc cx t^^ at later 

times. As a result, we expect a hardening spectrum at 
t < 2 days, and a softening spectrum at i > 2 days in 
the 100 TeV - 10 PeV range (Fig. 2). The peak time 
Thad is determined by equating E^"^ with the typical 
neutrino energy E^yP ~ 0.03^;^'^ [H]. We have Thad ~ 
2 X 10^ s V-if^33^'^f^sN^-iMl^l^i and the corresponding 



peak energy E^^^^^^ 



300 PeV 71^/^/"^/^/?^/" ... 
juu rev ?7_j Psn,-i^"sn,i 

which agree with Figs. 1 and 2. Note that cases of 

Tliad ^ 7em arc considered. The main contribution 

comes from the cosmic rays produced at t ^ Thad- Since 

El(j)^{t) (X E'^^{t) cx E^'' 1^, the neutrino fluence per 

flavor around the peak time is roughly estimated as ^ 

10-4 2 r)-2 r r 3/4„-7/4^3/4 -1/4 

iU tig cm i^g Mpc-/mcsjsup'/_i M33 A'SN.-l-^^SN,! ■ 

The total expected muon event rates (above 10°'^ TeV) 



-1/4 



by IceCube is ~ 2 D 



-2 

5 Mpc 



events in two days, which 



will be more than the atmospheric neutrino-induced 
event rates within 1°, Njf"^ ^ lO"^-^ events/day. Mag- 
netars arising at distances closer than 5 Mpc would yield 
higher fluxes observed as neutrino multiplets, which al- 
low us to recognize them as signals without coincident 
detections with photons and even to see the characteris- 
tic soft-hard-soft behavior. Since the magnetar birth rate 
is ~ 10^'^ yi'-^ galaxy-^, the probability to encounter a 
birth is non- negligible. From the number of local galax- 
ies, we expect ^ 0.02 — 0.05 yr"^ for the birth of magne- 
tars within 5 Mpc [lB| • 

One may expect an additional radiation field, leading 
to pry neutrinos in addition to pp neutrinos. For exam- 
ple, if the magnetar wind drives the SN explosion in its 
birth , a significant fraction of the outflow energy may 
be dissipated as radiation via the shocks. (The radia- 
tion field can also be expected in case of GRB jets in 
the star [l3|-) Therefore, we also show the case where 
the radiation field is included. In Figs. 1 and 2, the 



case for kT^ 



0.4 keVe\/'Elil,,f3, 



exp,5lA^SN.''-1^3 ^ ^^^^ 



shown. Here i^oxp is the outflow energy and is the 
radiation efficiency. When the radiation field exists in 
the SN ejecta, the previous expression of /mcs should 
be replaced with f^cs niin(l, max(/pp, /p^)), where 
the effective optical depth for the photomeson produc- 
tion, /p-y, is roughly estimated as fp^ w Kp^Cp^n^AsN — 
3/4 p3/4 ,q~5/4 ^-5/4 ^^.q^^j-^^ ^j^g A-resonancc en- 



OOU fc-y -C/oxp,5lA''SN,-l''4 

ergy of Ea ^ 2.4 PeV 



-l/4p-l/4 ^3/4 ,3/4 



0.2, cTp 



5 X 10 cm at the A-rcsonance. 



*^cxp,5i/?snViC- Here, 
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FIG. 2: Neutrino lightcurves corresponding to Fig. 1 at 
100 TeV (dashed line), 1 PeV (dotted line), 10 PeV (dotted- 
dashed line), 100 PeV (double dotted line), and the ratio of 10 
PeV fluence to 100 TeV fluence (thick line). Cooling of mesons 
and muons is important ai t < 2 days, while the amount of 
accelerated protons decreases with time. Thick/thin lines are 
for the cases without/with the radiation field. 



Correspondingly, the expression of /sup includes the cool- 
ing of mesons and muons due to interactions with pho- 
tons as well as their hadronic cooling. Following Ref. 
[T^ . neutrino spectra are numerically calculated, taking 
into account the radiation field. Although the radiation 
field can change spectra as a result of the difference in 
the meson multiplicity, we may expect that the total en- 
ergy fiuence around the peak energy and the qualitative 
feature are similar. 

Next, let us consider the sum of neutrinos from individ- 
ual magnetars, i.e., the cumulative neutrino background. 
The typical magnetar rate would be ~ 10 % of core- 
collapse (CC) SN rate, i?sN(0) - 1.2 x lO'^ Gpc"^ yr"! 
[U, . Possibly, the birth rate of fast rotating magnetars 
may be comparable to that of HNe that may be powered 
by magnetars, implying i?HN(0) ~ 2 x 10'^ Gpc"^ yr"^ 
[l8|. By using our numerical results, the cumulated fluxes 
can be estimated as [13, [3 

- 3 X 10-^ GeVcm-2 str" V„,es/sup?/!i/i^3^ 



„ of /\r^3 /gco /z ^mag(O) 

^ /-'SN,-1-'"SN,1 



0.5 3 1.2 X 10^ Gpc"^ 



--(4) 



where fg^o is the fraction of the magnetars with the pre- 
ferred geometry for ion acceleration Q and fz expresses 
the contribution from the high redshift sources [ij, [l^ . 

The numerical results are shown in Fig. 3, where the 
birth rate evolution is included with the SFR2 model 



for magnetars [l4l |. The muon event rates are N^j, 
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events/yr for i?niag = 0.1i?sN and A^^ ~ 4 events/yr for 
-Rmag = -Rhn, respectively. Note that cross-correlation 
studies between neutrinos and CC SNe/HNe can be 
important, but we may expect one time- and space- 
coincident event among ~ 10'*-^ magnetar births. How- 
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FIG. 3: The cumulative (//g + + v^^ + v^^ + u-r + Vt) back- 
ground from newly born magnetars. Thick solid line: Mag- 
netar (10% of CC SNe) shows the case where 10% of CC 
SNe bear fast rotating magnetars, with Msn ~ IOM0. Thin 
solid line: Magnetar (Hypernovae) is a case where the fast 
rotating magnetar rate is comparable to the HN rate, with 
Msn ~ Mq. For comparison, various GRB neutrino back- 
grounds for the GRB-UHECR hypothesis are also shown. 
Thick dashed line: HL GRB, prompt neutrinos from HL 
GRBs [11, [il. Thin dashed line: LL GRB, prompt neutri- 
nos from LL GRBs [13]. WB: The Waxman-Bahcall bounds 
are shown as benchmarks [3. Atm: the conventional atmo- 
spheric neutrino background. A ACDM cosmology {Qm ~ 
0.3, flA ~ 0.7; Ho — 71 kms^'^ Mpc~^) and Zmax ~ 11 is used. 



ever, from Fig. 3, we still expect future constraints on 
the "diffuse" neutrinos arc important to test the mag- 
netar scenario since expected neutrinos have rather high 
energies. 

We consider the possible beaming effects of the cosmic 
rays, SNR puncture and disruption by the winds, since 
these can affect the escape and detectability of the cosmic 
rays. Cosmic rays themselves may be beamed, enhanc- 
ing the neutrino signals from individual sources by the 
beaming factor while the background is unchanged. In 
addition, if the PNS winds are significantly collimated, 
they may puncture the stellar envelope, leading to long 
GRB jets Only cosmic rays that are emitted along 
the penetrating jets can escape without depletion. How- 
ever, even if cosmic rays are beamed along the jet, we still 
expect neutrino production when jets are choked rather 
than successful |17| . 

The disruption resulting in the formation of supershcUs 
expanding into the interstellar medium was discussed in 
Ref. [^. Such phenomena have never been observed in 
CC SNe/HNe, but we discuss them for the sake of gen- 
erality. The effect of the disruption by energy shedding 
may be characterized by a dumpiness factor, C = Sp/p. 
The no shedding case corresponds to C = 1. In clumpy 
SNRs, neutrinos are produced when clumps lie along the 
line of sight. While values of C are uncertain, as an 
example C ~ S'^ implies that the probability to see neu- 



trinos is ~ C ^ 1/9, and the background is similarly 
reduced. 



IV. IMPLICATIONS AND DISCUSSIONS 

The birth of fast rotating magnetars is a common sce- 
nario discussed in connection with HNe and GRBs j^, 
although the origin of their strong magnetic fields is con- 
troversial (e.g., dynamo vs fossil fields). Here we have 
shown that high-energy neutrino signals can serve as the 
smoking gun signal announcing the birth of fast rotating 
magnetars. A suppression of the highest-energy neutri- 
nos at t < Thad ~ a few days implies that cosmic-ray 
acceleration occurs inside the SNR, and a characteristic 
soft-hard-soft behavior is expected. Although our pre- 
dicted fluxes are below the current observational limits, 
they can be tested by IccCubc, KM3Nct, ARIANNA and 
Auger in the near future. Even non-detections would pro- 
vide useful constraints on this magnetar scenario. 

There arc also other possible ways to produce neutri- 
nos. First, at earlier times (< days), there may be ra- 
diation fields with the temperatures of 0.1 — 1 kcV (at 
rsN ~ 10^'^'^ cm), in the shocked stellar ejecta ahead 
of the wind [H, [l^ or in the cavity The former 

case is also demonstrated in this work. Collimated winds 
launched at very early times (< lO'^ s) may become suc- 
cessful jets such as GRB jets. Then, a fraction of the 
cosmic rays interact with late internally dissipated or 
external-shock photons, making other very high-energy 
neutrino signals, similarly to the case of GRBs [l^ . 

Cosmic rays and neutrinos could be expected possibly 
also in normal pulsars where Thad ^ ?em is anticipated. 
When they have weaker magnetic fields but rapid rota- 
tion speeds, the peak time of the energy fluence can be 
later than days, and then neutrino emission lasts longer. 
Since it was not the focus of this work, see Rcfs. for 
comprehensive discussions. 

Magnetar neutrinos may be useful for revealing the 
possible connection between magnetars and UHECRs. 
Since our purpose here is not an explanation of UHECRs, 
wc make only two brief comments on this possibility: (a) 
only a fraction of magnetars can be UHECR sources, 
and (b) the Auger spectrum seems to conflict with Eq. 
([^ |2l| . Concerning point (a) , one may think of many 
possible reasons such as initial rotation rate, geometry, 
puncture or disruption of the SNR. Interestingly, the rate 
required in the magnetar scenario is comparable to the 
HN rate [2lj . Point (b) requires more careful consid- 
eration, but we may expect a realistic spectrum to differ 
from Eq. depending on the detailed mechanism. For 
example, if test particles are stochastically accelerated 
by waves, dN/dE oc E~'^ can be expected [l3|. Although 
this requires further investigations beyond the scope of 
this work, the use of Eq. ^ is sufficient for demon- 
strations. Other cases can easily be predicted, once the 
cosmic-ray spectrum is given. 



5 



Acknowledgments 

K.M. thanks K. Kohri for helpful discussions. We ac- 
knowledge support by a Grant-in- Aid from JSPS, and by 
a Grant-in- Aid for the Global COE Program "The Next 



Generation of Physics. Spun from Universality and Emer- 
gence" from MEXT (K.M.), by NASA NNX08AL40G, 
NSF PHY -0757155 (P.M.) and NASA NNG05GB67G 
(B.Z.). 



[1] P.M. Woods and C. Thompson, Compact Stellar X-ray 
Sources, eds. W.H.G. Lewin and M. van der Klis, Cam- 
bridge Astrophysics Series, 39, 547 (2006); A.K. Hard- 
ing and D. Lai, Rep. Prog. Phys. 69, 2631 (2006); S. 
Mereghetti, Astron. Astrophys. Rev. 15, 225 (2008). 

[2] R.C. Duncun and C. Thompson, Astrophys. J. 392, L9 
(1992); C. Thompson and R.C. Duncan, Astrophys. J. 
408, 194 (1993). 

[3] P. Meszaros, Rep. Prog. Phys. 69, 2259 (2006); B. Zhang, 
Chin. J. Astron. Astrophys. 7, 1 (2007). 

[4] T.A. Thompson, P. Chang, and E. Quataert, Astrophys. 
J. 611, 380 (2004); N. Bucciantini et al, Mon. Not. 
R. Astron. Soc. 368, 1717 (2006); B.D. Metzner, T.A. 
Thompson, and E. Quataert, Astrophys. J. 659, 561 

(2007) ; S.S. Komissarov and M.V. Barkov, Mon. Not. 
R. Astron. Soc. 382, 1029 (2007). 

[5] J. Arons, Astrophys. J. 589, 871 (2003). 

[6] G. Ghisellini et al., Mon. Not. R. Astron. Soc. 390, L88 

(2008) . 

[7] W. Bednarek and R.J. Protheroe, Phys. Rev. Lett. 79, 
2616 (1997); R.J. Protheroe, W. Bednarek, and Q. Luo, 
Astropart. Phys. 9, 1 (1998); W. Bednarek, Mon. Not. 
R. Astron. Soc. 331, 483 (2002); J.H. Beall and W. 
Bednarek, Astrophys. J. 569, 343 (2002). 

[8] B. Zhang et al, Astrophys. J. 595, 396 (2003); T. Herpay 
et al, J. Cosmol. Astropart. Phys. 08, 025 (2008). 

[9] J. Ahrens et al, Astropart. Phys., 20, 507 (2004); 
U.F. Katz, Nucl. lustrum. Methods Phys. Res., 
Sect. A 567, 457 (2006); V. Van Elewyck et al., 
arXiv :astro-ph/06 1 273 1 
[10] P. Chen, T. Tajima, and Y. Takahashi, Phys. Rev. Lett. 
89, 161101 (2002); M. Hoshino, Astrophys. J. 672, 940 



(2008); Y. Kuramitsu et ai, Astrophys. J. 682, L113 

(2008) . 

[11] Y.A. Gallant and J. Arons, Astrophys. J. 435, 230 
(1994). 

[12] P. Blasi, R.I. Epstein, and A. Olinto, Astrophys. J. 533, 
L123 (2000). 

[13] R.S. Fletcher, T.K. Gaisser, P. Lipari, and T. Stanev, 

Phys. Rev. D 50, 5710 (1994); S.R. Kelner, F.A. Aharo- 

nian, and V.V. Bugayov, Phys. Rev. D 74, 034018 (2006). 
[14] K. Murase, Phys. Rev. D 76, 123001 (2007); Phys. Rev. 

D 78, 101302(R) (2008). 
[15] S. Razzaque, P. Meszaros, and E. Waxman, Phys. Rev. 

D 68, 083001 (2003); S. Ando and J.F. Beacom, Phys. 

Rev. Lett. 95, 061103 (2005). 
[16] S. Ando, J.F. Beacom, and H. Yuksel, Phys. Rev. Lett. 

95, 171101 (2005). 
[17] P. Meszaros and E. Waxman, Phys. Rev. Lett. 87, 171102 

(2001). 

[18] P. Madau, M. Delia Valle, and N. Panagia, Mon. Not. R. 

Astron. Soc. 297, L17 (1998); D. Guetta and M. Delia 

Valle, Astrophys. J. 657, L73 (2007). 
[19] E. Waxman and J. BahcaU, Phys. Rev. Lett. 78, 2292 

(1997); Phys. Rev. D 59, 023002 (1998); K. Murase and 

S. Nagataki, Phys. Rev. D 73, 063002 (2006); Phys. Rev. 

Lett. 97, 051101 (2006). 
[20] S.E. Woosley, N. Langer, and T.A. Weaver, Astrophys. 

J. 411, 823 (1993). 
[21] J. Abraham et al, Phys. Rev. Lett. 101, 061101 (2008); 

K. Murase and H. Takami, Astrophys. J. 690, L14 

(2009) . 



